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ABSTRACT 


Large  deformation  elastic-plastic  buckling  loads  are  ob- 
tained for  axisymmetric  spherical  caps  with  initial 
imperfections.  The  problem  formulation  is  based  on  equilib- 
rium equations  in  which  the  plastic  deformation  is  taken  as 
an  effective  plastic  load.  Both  perfectly  plastic  and  strain 
hardening  behavior  are  considered.  Strain  hardening 
is  represented  by  the  Prager-Ziegler  kinematic  hardening 
theory,  so  that  the  Bauschinger  effect  is  accounted  for. 
Solutions  of  elastic-plastic  circular  plates  and  spherical 
caps  are  in  good  agreement  with  previous  results.  For  the 
spherical  cap  it  was  determined  that  both  initial  imperfec- 
tion and  plastic  deformation  have  the  same  effect  of  reducing 
buckling  capacity;  as  the  magnitude  of  the  imperfection  in- 
creases, the  influence  of  plastic  deformation  becomes  less 
important.  It  is  also  found  that  the  geometric  parameter  A, 
which  is  used  as  an  important  factor  in  elastic  response,  be- 
comes meaningless  for  the  elastic-plastic  buckling  analysis 
of  spherical  caps. 
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INTRODUCTION 


* 


To  realistically  predict  buckling  loads  for  spherical 
caps,  consideration  should  be  given  to  the  effect  of  geome- 
tric and  material  nonlinearities  as  well  as  initial 
imperfections,  since  all  three  factors  play  significant  roles 
in  reducing  spherical  cap  buckling  capacity. 

The  axisymmetric  buckling  of  clamped  shallow  spherical 
shells  were  studied  by  Huang  [1]  and  Weinitschke  [2],  To  re- 
duce the  gap  between  experiment  and  theory,  they  introduced 
asymmetrical  deflections  along  the  circumferential  direction 
and  obtained  buckling  pressures  for  different  cap  geometric 
parameters  by  solving  eigenvalue  problems. 

Kho jasteh-Bakht  [3]  analyzed  small  deflection  elastic- 
plastic  shells  of  revolution  by  the  finite  element  method. 
Marcal  extended  the  work  of  [3]  to  include  large  deforma- 
tion terms  [4].  An  isotropic  strain-hardening  is  assumed  in 
these  two  references  and  it  is  found  in  Ref.  [4]  that  the 
plastic  yielding  has  the  significant  effect  of  reducing  the 
buckling  pressure  of  shell  structures.  Various  solution  pro- 
cedures for  the  large  deformation  elastic-plastic  analysis  of 
shell  structures  were  discussed  in  Ref.  [5].  Shells  of  revo- 
lution under  cyclic  loading  were  examined  in  Ref.  [6] . 
Bushnell  [7]  employed  a variational  finite-difference  ap- 
proach to  investigate  the  influence  of  plastic  deformation  on 
a shell  structure. 

Initial  imperfections  in  the  shell  shape  sometimes  ac- 
counted for  the  discrepancy  between  experimentally  observed 
buckling  loads  and  theoretically  calculated  values  [8]. 
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Hutchinson  [9]  utilized  the  Koiter  approach  to  determine  the 
approximate  asymmetric  buckling  load  of  spherical  shells  with 
initial  imperfections.  Asymmetric  buckling  loads  for  spheri- 
cal caps  with  asymmetric  imperfections  were  obtained  in  Refs. 
[10,11]  by  directly  solving  two  dimensional  governing  differ- 
ential equations. 

It  is  the  purpose  of  this  work  to  study  the  effect  of 
initial  imperfections,  and  geometrical  and  material  nonlinear- 
ities on  the  axisymmetric  buckling  of  spherical  caps.  It  is 
hoped  that  the  result  obtained  may  provide  a better  under- 
standing on  the  axisymmetric  buckling  nature  of  spherical 
caps.  The  analysis  of  the  proposed  title  problem  can  be  pro- 
ceeded by  incorporating  the  appropriate  plasticity  theory 
into  an  existing  computer  code  [12], 

Selection  of  a plasticity  theory  is  very  important.  The 
theory  should  be  simple  in  mathematical  expression,  easy  to 
apply  and  capable  of  representing  the  material  property  such  as 
the  Bauschinger  effect.  In  accordance  with  these  criteria, 
the  von  Mises  initial  yield  condition,  the  incremental  flow 
theory  [13]  and  the  Prager-Ziegler  kinematic  hardening  theory 
[14,15]  are  used  in  this  study. 

The  problem  formulation  is  based  on  equilibrium  equa- 
tions, treating  the  plastic  deformation  as  effective  plastic 
loads  which  are  lumped  together  with  actual  applied  loads 
[16].  Governing  equations  are  replaced  by  finite  differences 
and  resulting  equations  are  then  solved  by  the  nonlinear  re- 
laxation method  [17]. 

In  the  next  section,  the  governing  equations  based  on 
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equilibrium  considerations  are  given.  This  is  followed  by 
the  description  of  plasticity  relations  and  the  general  solu- 

» 

tion  procedure.  To  test  the  validity  of  the  plasticity  theory 
adopted,  example  problems  are  solved  and  comparisons  with 

k 

other  solutions  are  made.  Large  deformation  elastic-plastic 
buckling  loads  are  then  obtained  for  spherical  caps  with  in- 
itial imperfections.  Conclusions  and  a general  discussion 
are  given  in  the  final  section. 

* 
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PROBLEM  FORMULATION 

The  geometry  of  a spherical  cap  is  shown  in  Fig.  la,  in 
which  H is  the  central  height  and  R the  shell  radius;  a is  the 
base  radius;  W(r)  and  U(r)  are  the  displacement  components 
along  normal  and  tangential  directions,  respectively,  and 
Wi(r)  is  the  initial  imperfection;  q is  the  applied  uniform 
pressure.  The  undeformed  shape  of  the  perfect  shell  can  be 
adequately  described  by 

Z = H[l-(r /a)2]  (la) 

and  the  radius  of  curvature  of  the  shell  is  approximated  by 

R = a2/ 2H  (lb) 

where  r is  the  radial  coordinate. 

Fig.  lb  shows  the  membrane  forces  N and  NQ , the  trans- 

r o 

verse  shear  Q and  the  moments  M and  M. . Equilibrium  of 
r r d 

moments  requires 


(rMr) ’-MQ-rQr  = 0 (2) 

And  equilibrium  of  stress  resultants  along  radial  and  normal 
directions  provides 


(rNr) ' -N0  = 0 (3) 

lrNr(Wf-Z) ’+rQr) ’+rq  = 0 (4) 

where  ( )'=3(  )/9r  and  W^=W+W^ . Note  that  the  nonlinearity 
has  been  introduced  in  Eq.  (4)  by  considering  the  influence 
of  W^. 

Elimination  of  Qr  in  Eqs . (2)  and  (4)  yields 


- M'+N  (Wl'+ 
r 0 r f 


1 Wf 

r)+Vf- 


= 


(5) 
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* 


Eqs . (3)  and  (5)  are  the  basic  equations  for  the  analysis  of 


axisymmetric  spherical  caps, 


Stress-Strain  and  Strain-Displacement  Relations 


For  a shell  deforming  into  plastic  range,  the  strain  in 


a point  within  the  thickness  can  be  expressed  by 


fe}  = {e  }+{e^*} 


where  {e},  {ee}  and  {e^3}  are  the  total,  elastic  and  plastic 


strain  vectors,  respectively. 


Furthermore,  the  total  strain  can  be  considered  as  the 


sum  of  the  membrane  and  bending  components: 


{e}  = {e}  + z{i<} 


where  z is  the  vertical  coordinate  through  the  shell  thick- 


ness (Fig.  lb) . The  membrane  and  bending  strains  are  related 


to  displacements  by 


e * U'-  ^ + i(W')2+W’W? 

i K Z 1 


„ U w 

e0  r " R 


< = — W" 

r 


kq  = - — 
0 r 


Note  that  the  elastic  components  of  strains  are  the  only 


strains  which  can  be  related  to  stresses  by  Hook's  law: 


(o)  = [E] ( (e}-{ep) ) 


where  [E],  the  elastic  strain  to  stress  transformation  matrix, 


is  given  as 


1 

1-v  Lv  lJ 
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in  which  E is  Young's  modulus  and  v is  Poisson's  ratio. 

Membrane  stress  resultants  and  bending  moments  are  ob- 
tained by 


{N} 

{M} 


{a}dz 


{a} zdz 


(ID 

(12) 


Substituting  Eqs.  (6-9)  into  Eqs.  (11)  and  (12),  we  obtain 
the  membrane  forces 


where  the  effective  plastic  membrane  forces  are 


Np 

r 

; = [E] 

r h/2 

(eP 

1 r 

nP 

0 

j « 

<N 

\ 

r* 

1 

ep 
»e  0 

and  the  moments 


(13) 


(14) 


n ■ »r  i fri  - iM? 

K>  L lJ  W K 


(15) 


3 2 

where  D=Eh  /12(l-v  ),  and  effective  plastic  moments  are 


MP 

r 

= [E] 

r h/2 

eP 

r 

mP 

0 

■'-h/2 

ep 

e0 

zdz 


(16) 


Governing  Equations 


In  terms  of  displacements,  the  governing  equation  in- 
volving the  major  displacement  U is  obtained  by  substituting 
Eqs.  (13)  and  (14)  into  Eq.  (3)  : 


U"  + 


ir 

r 


- \ +G(W) 
r 


(17) 


where 


rr 
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1-v 


G (W)  = F ' (W) +vF ' (W) + 

r 0 r 


(Fr-Fe) 


Fr(W)  = -W/R+(W' )2/2+W'W| 
F0(W)  = -W/R 

F' (W)  = -W,/R+W,W"+W,WV+W"W! 
r 11 


(18) 


F 0 ( w ) = -W'/R 

and  , an  effective  plastic  load,  is  expressed  in  terms  of 
effective  plastic  membrane  forces  by 

qP  = (N^)  ' +N^/ r-N0/ r (19) 

The  governing  equation  involving  the  major  displacement 
W is  also  obtained  by  substituting  Eqs.  (13-16)  into  Eq . (5) : 


4 Fh 

DV  w-  -^(er+ve0)  (W£+l/R)-  JVoQ.  _r,  v..f/ 


Eh 


( e +ve  ) (W'/r+l/R) 


1-v 


P_„P 


= q-q^-q 


(20) 


4 2 2 2 

where  V =V  (V  ) and  V ( )=(  )"+(  ) ' /r;  the  membrane  strains 
Lr'e0  are  defined  in  Eqs.  (8),  effective  plastic  loads  qP  and 
qP  are  given  as 


(21) 


qP  = N P(W£+1/R. ) +nP  (W^/r+l/R) 

qP  = (mP)  ”+2(mP)  7r-(M0p)  • /r  (22) 

Eqs.  (17)  and  (20)  are  two  fundamental  governing  equa- 
tions in  terms  of  displacements  for  the  present  analysis. 

Boundary  Conditions 

At  shell  apex,  the  nature  of  axi-symmetry  requires  that 

W' (0)  = 0 (23) 


U(0) 


0 


(24) 
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Along  the  outer  edge  (r=a) , if  the  cap  is  clamped: 

U ( a)  = W (a)  = W'  (a)  = 0 (25) 

On  the  other  hand,  if  the  cap  is  simply  supported,  it  re- 
quires that 

U ( a)  = W (a)  = 0 (26a) 

and  that  Mr(a)  in  Eq.  (15)  to  be  zero,  i.e. 


dw  + V dw 
dr2  r dr 


-M^(a)f  r =a 


(26b) 


where  is  defined  in  Eq.  '16). 


Nondimensional  Forms 


For  convenience,  the  following  nondimensional  quantities 
are  introduced 


x = r/a 

2 2 2 
A = m a /Rh 

( ) ' = 3(  )/3x 

u = aU/h2 
w = W/h 


4 2 

ni  = 12  (1-v  ) 

2 2 2 

qcr  = 4Eh  /R  in 
p = (2  7) 

wi  = wi/h 


where  qcr  is  the  classical  buckling  pressure  of  a complete 
spherical  shell  of  the  same  radius  of  curvature  and  thickness 
With  the  adoption  of  Eqs.  (27),  the  nondimensional  forms  of 
Eqs.  (17)  and  (20)  become 


u"  + 


u' 


u j.  i,  \ (1-v2)  a ^ p 

j +g(w)  = L qt' 

x Eh  i 


(28) 


and 
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vVi2(cr+vr0)(w»+  4)-i2<Vver,(“i  ♦ 4) 

' m / \ m / 


,4  4 4 

= 4 T-  P-  ^(qf+q?) 

m Eh  ^ J 


(29) 


where  g (w)  and  f (w)  terms  are  nondimensional  counterparts  of 
similar  terms  in  Eq.  (18): 

g(w)  = f^(w)+vf • (w)+(l-v) [fr(w)-f  (w)]/x 
, 2 


f (w) 
r 


f0(w) 


f ' (w) 
r 


— ^ w+  ^(w')2+w'w^ 


A2 

T w 

m 


(30) 


A2 

7 


w'  +w' w"+w'  w'.'+w"w! 


f Q (w)  = - — w' 


m 


an<3  £r»  fcg  are  nondimensional  quantities  of  membrane  strains 
er  and  eQ  in  Eq.  (8) 


T 

m 


w+  i(w' ) 2+w,w'.' 

2 l 


(31) 


- _ u A 

C n — — “ — W 

Ox  2 
m 


W 
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PLASTICITY  RELATIONS 

As  indicated  by  Ziegler  [15] , the  behavior  of  a strain- 
hardening material  can  be  described  by 

(a)  an  initial  yield  condition,  defining  the  elastic 
limit  of  the  material, 

(b)  a flow  rule,  relating  the  plastic  strain  increment 
with  the  stress  and  the  stress  increment, 

(c)  a hardening  rule,  specifying  the  subsequent  yield 
condition  in  the  course  of  plastic  flow. 

Because  of  complexity  associated  with  plastic  deforma- 
tion, the  laws  governing  the  behavior  of  materials  in  the 
plastic  range  have  not,  as  yet,  reached  a level  of  general 
acceptance.  A judicious  choice,  among  all  available  plasti- 
city theories,  should  consider  both  the  simplicity  in  mathe- 
matical expression  and  the  proper  representation  of 
experimentally  observed  material  behavior. 

Based  on  these  considerations,  the  von  Mises  yield  func- 
tion, which  describes  a smooth  surface  in  stress  space  and 
represents  a simple  mathematical  function,  is  chosen  as  the 
initial  yield  condition.  The  flow  rule  of  von  Mises  and  the 
Ziegler-Prager  kinematic  hardening  rule  are  also  selected  in 
this  study.  It  is  important  to  note  that  this  hardening  rule 

I 

does  account  for  the  Bauschinger  effect. 

In  a 9-space  stress  field  with  origin  0 (Fig.  2) , the 
von  Mises  initial  yield  surface  can  be  described  by 

F(°ij)  = k2  = constant  (32) 


For  an  initially  isotropic  material,  the  form  of  the  function 


F is  invariant  with  respect  to  a rotation  of  the  stress  state. 

The  Prager's  hardening  rule  [14]  assumes  that  during 
plastic  deformation,  the  yield  surface  moves  in  a translation 
without  changing  its  shape  and  hence  the  subsequent  yield  sur- 
face takes  the  form 

F(aij"Ctij)  = k2  (33) 

where  represents  the  total  translation  of  the  yield  surface 
center  which  is  a measure  of  the  degree  of  work  hardening.  In 
the  space  o^,  ou  ^ is  the  position  vector  of  the  yield  surface 
center  C which  before  plastic  deformation  takes  place  is  lo- 
cated at  the  origin  (Fig.  2) . 

The  flow  rule  of  von  Mises  gives  the  following  expression: 

de? . = 7t— — — dA  , dA>0  (34) 

i n do  . . 

J ID 

which  indicates  that  the  plastic  strain  increment  de? . lies 
in  the  qxterior  normal  of  the  yield  surface  (33). 

The  Ziegler's  modification  [15]  of  Prager's  hardening 
rule  suggests  that  the  surface  (33)  moves  in  the  direction  of 
the  radius  connecting  its  center  with  the  stress  point  (Fig.  2) 
dai j = (o^-a^du  , dy>0  (35) 

dA  in  Eq.  (34)  and  dy  in  Eq.  (35)  are  to  be  determined. 

Determination  of  dy  is  based  on  the  condition  that  stress 
point  always  remains  on  the  yield  surface  in  plastic  flow. 

This  condition,  in  fact,  states  that  for  an  infinitesimal  in- 
crement of  loading  the  vector  do^ j-d°u  j must  be  orthogonc ’ to 
the  outer  normal  to  the  yield  surface: 
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Substituting  (35)  into  (36)  yields 


dy  = 


(SF/3<Wd°k<t 

7 a -a  ) 9F/9o 
mn  mn  mn 


(37a) 


and  hence 


d0tij  To  -cT*)  9F/9a  (aij~aij) 

J mn  mn  ' mn  J J 


OF/30kl)aokt 


(37b) 


mn  mn 


mn 


According  to  Refs.  [14,15],  9 A in  (34)  can  be  obtained 


by  assuming  da^=cde:\  and  putting  this  relation  into  (36): 
dX  = I (9F/8ak&)dok& 

c (3F/3a'  ) (9F/9a  ) (38a) 

mn  mn 

and  hence 


_ 1 <3F/3ok0)dOk«.  3F 

ij  * c (3F/3ainnH3F/3(,mn)  3o.  . 


(38b) 


where  c is  the  hardening  coefficient  and  can  be  determined 
from  the  uniaxial  stress-strain  relation. 

Now,  let  us  write  Eq.  (33)  in  terms  of  3 principal 
stresses : 

2 


f = F-k 


= ^-[  (o1-a2)2+(a2-o3)2+(o3-a1)  2]-o3  = 0 


(39) 


where  is  the  yield  stress  in  uniaxial  tension  and  = 

°l-al'  a 2~a 2~a 2 ' a 3=a 3_a 3 • 

For  the  case  of  plane  stress  (Fig.  3) , we  have 


°3  = “3  = °3  = 0 


and  Eq.  (39)  reduces  to 


—2  — — —2  7 

f = o .-o.o~+oZ-o  L = 0 
1 1 2 2 y 


(40) 


(41) 


By  putting  Eq.  (39)  into  Eq.  (38b)  and  using  relations  (40-41), 


■mmmh MM 


From  Hook's  law,  we  have 


I40l] 

E 

'1  v“ 

Ael 

|/So2| 

i 2 

1— V 

.V  1. 

Ae2 

It  is  also  noted  that 


where 


'11 


|fle! 

AS1 

K 

'Ae2 

Ae2 

K 

(42) 

into  (44) 

provides 

fA°l 

rCll 

C12  "I  jAe 

> 

Q 

to 

-C21 

C22J  (Ae 

(D+E 

s2) 

C12  = i!(Dv-ESlS2)  = C21 


(42) 


| c (43) 


(44) 


(45) 


(46) 


C22  = §<D+E  sj) 

U = D(1-v2)+E(S^+2vS1S2+S2) 

For  a given  displacement  field,  {Ae}  are  obtained  from  incre- 
mental forms  of  Eqs . (7)  and  (8).  Eqs . (46)  are  the 

stress-strain  relations  during  the  course  of  plastic  flow 
(loading).  Otherwise,  Eqs.  (44)  should  be  used  for  all  stress 
computations . 
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Up  to  this  point,  the  hardening  coefficient  c in  Eq. 

(43)  is  the  only  constant  yet  to  be  decided.  If  the  structure 
is  in  a state  of  uniaxial  stress,  the  stress-strain  relation 
is  the  same  as  that  obtained  from  tension  or  compression  tests. 
The  expression  for  the  hardening  coefficient  c can  be  readily 
obtained  from  Eq.  (42)  by  setting  AeP=AeP,  Aa=Aa1: 


or 


1 

c 


D = 


3 Ae^ 
2 Aa 


Aa 

Ae1 


(47a) 


(47b) 


This  equation  shows  that  the  value  of  D is  equal  to  the  slope 
of  the  uniaxial  stress-plastic  strain  curve. 

From  Eqs.  (47) , the  values  of  the  hardening  coefficient 
for  two  special  cases  may  be  specified  here:  (i)  for  an 
elastic-ideally  plastic  material,  c=D=0,  (ii)  for  a linear 
hardening  material,  D=EEt/ (E-Et) , where  Efc  is  the  tangent 
modulus.  Uniaxial  stress-strain  curves  for  these  two  types 
of  hardening  material  are  given  in  Fig.  4. 

Nonlinear  Hardening 

Also  given  in  Fig.  4 is  the  stress-strain  curve  for  a 
nonlinear  hardening  material.  Its  D value,  instead  of  being 
a constant,  depends  on  the  state  of  stress.  One  way  of  deal- 
ing with  this  rather  complicate  situation  is  based  on  the 
Ramberg-Osgood  representation  of  a uniaxial  stress-strain 
curve  [18]: 

n-1 


o . 3o 

e = E 7E 


0.7 


(48) 


where  n = 1+ 


log (17/7) 


log(°0.7/°0.85) 

e is  the  total  strain,  E is  the  slope  of  the  linear  portion 
of  the  stress-strain  curve,  and  oQ  7 and  aQ  g5  are  the  stress- 
es at  which  the  curve  has  secant  moduli  of  0 . 7E  and  0.85E, 
respectively . 

It  is  understood  that  the  nonlinear  term  in  Eq.  (48)  is 
the  plastic  strain.  D value  of  this  material  can  be  obtained 
as 

n-1 


D » ¥ 
3n 


0.7 


(49a) 


We  may  generalize  this  equation  to  a multiaxial  state  of 
stress  [19,20] : 

o. 


D = 


7E 

3n 


0.7 


n-1 


(49b) 


where  a,  the  effective  stress,  is  defined  as 


a = J~\ r a ! . o ! . 

2 i]  i] 


(50a) 


in  which  ^ is  the  deviatoric  stress.  For  the  case  of  plane 
stress , 


2 2 
a — a i~° 2+a  2 


(50b) 


Loading  Criteria 

In  the  incremental  solution  procedure  for  elastic-plastic 
problems,  in  addition  to  the  constitutive  relations,  it  is 
necessary  to  have  a loading/ unloading  criterion.  For  this 
purpose,  let's  define  f=Of/3oi^)  do^,  where  f has  been  de- 
fined in  Eq.  (39).  Loading,  unloading  and  neutral  loading 
are  associated  with  the  plastic  state  f=0,  and  are 
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characterized  by  f>0,  f<0  and  f=0,  respectively.  When  load- 
ing or  neutral  loading  takes  place,  Eq.  (46)  must  be  applied. 

On  the  other  hand,  Eq.  (44)  must  be  used  when  unloading  occurs. 


SOLUTION  PROCEDURE 


For  convenience,  a simple  flow  chart  is  sketched  in  Fig. 

5 to  explain  the  general  solution  procedure.  The  entire  pro- 
cess is  divided  into  two  major  loops,  namely,  the  elastic 
solution  and  material  property  loops. 

In  the  elastic  solution  loop,  all  material  properties  are 
held  constant.  Consequently,  the  effective  plastic  loads, 
q?  in  the  two  governing  equations  (28)  and  (29)  are  fixed  and 
combined  with  the  actual  externally  applied  load  q,  the  prob- 
lem is  thus  reduced  to  an  elastic  large  deformation  problem. 

Central  finite  differences  are  used  to  convert  Eqs.  (28) 
and  (29)  to  a discrete  system  of  equations,  and  the  nonlinear 
relaxation  technique  is  employed  to  solve  these  nonlinear 
equations.  For  details  of  solving  these  equations  by  the  non- 
linear relaxation  technique,  readers  are  referred  to  Ref. 

[17].  The  iteration  in  this  loop  is  considered  converged 
when  the  average  absolute  change  of  u and  w displacements  at 
all  points  is  less  than  0.0001. 

With  the  new  displacement  field  {u}  , the  material  prop- 
erties must  be  updated  so  that  the  (nonlinear)  stress-strain 
relation  can  be  satisfied  at  all  points  over  the  shell  surface 
and  through  the  thickness.  In  the  material  property  loop, 
the  loading  criterion  is  first  checked.  If  the  material  point 
is  in  an  unloading  situation  or  still  in  the  elastic  range. 
Hook's  law  is  used  and  the  computation  is  very  straight-for- 
ward . 

If  the  material  point  is  in  a loading  situation,  incre- 
mental strains  {Ae}  are  computed  from  Eq.  (7)  by  using 
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3c  k—1 

{Au}  = {u}  - (u)  (where  k is  the  number  of  the  material  property 

loop,  see  Fig.  5),  and  (Ao)  are  computed  from  Eq.  (46).  Having  obtained 
{Aa}  , the  new  values  of  {a},  {eG}  and  {eP}  can  readily  be  calculated. 
Effective  plastic  loads  q?  are  then  evaluated  from  integration  formulae 
(19),  (21)  and  (22)  through  the  use  of  Simpson's  rule  (9  thickness  points 
are  used  for  the  entire  computation  in  this  paper). 

The  material  property  is  considered  to  be  updated  if  average  absolute 
change  of  displacements  between  the  present  and  the  previous  material 
property  loop  is  less  than  0.001.  Otherwise,  the  iteration  goes  to  the 
elastic  solution  loop  again  and  the  entire  operation  is  repeated  until 
the  material  property  and  equilibrium  equations  are  satisfied. 


COMPARISON  WITH  EXISTING  SOLUTIONS 


To  test  the  validity  and  accuracy  of  plasticity  theory 
utilized  in  this  study,  three  small  deformation  circular 
plates  and  two  large  displacement  spherical  caps  of  various 
material  properties  are  chosen  as  example  problems. 

The  first  example  is  a uniformly-loaded  simply  supported 
circular  plate  (Fig.  6a) . Elastic-ideally  plastic  material 
is  assumed  with  the  value  of  the  fully  plastic  moment,  Mq, 
equal  to  4000  lb  in/in.  11  and  9 points  are  selected  along 
the  plate  surface  and  thickness,  respectively.  After  the 
elastic  limit  has  reached  the  load  was  increased  in  incre- 
ments of  4 to  5 psi. 

Results  are  displayed  in  Fig.  6a-6c;  also  shown  in  these 
figures  are  solutions  of  Ref.  [21].  Ref.  [21]  adopted  a fi- 
nite element  formulation  and  the  isotropic  hardening  rule. 
Despite  differences  in  the  method  of  solution,  these  two  sets 
of  results  are  virtually  identical. 

Also  shown  in  these  figures  are  solutions  from  the  elas- 
tic analysis.  A redistribution  of  displacements  and  moments 
as  a result  of  the  plastic  deformation  is  quite  apparent  in 
these  figures. 

The  same  example  problem  is  considered  again  with  the 
exception  that  the  material  property  is  of  linear  hardening 
behavior  as  Riven  in  Fig.  7a  and  7b.  The  same  surface  and  thickness  points 
arc  used  and  results  are  displayed  in  Fig.  7a  and  7b.  Also  plotted  in 
these  figures  are  elastic  solutions  and  those  calculated  by 
Kho jasteh-Bakht  (3).  He  used  a finite  element  analysis  and 
an  isotropic  hardening  rule.  Comparison  of  present  results 
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and  those  of  [3]  is  excellent.  Again,  the  effect  of  the  plas- 
tic deformation  on  displacements  and  moments  is  obvious  from 
these  figures. 

The  third  example  problem  is  a uniformly  loaded  clamped 
circular  plate  of  nonlinear  hardening  material  (Fig.  8a) . 

This  problem  was  originally  considered  in  Ref.  [22].  The  ma- 
terial stress-strain  curve  given  in  Ref.  [22]  can  be  described 

7 

by  the  following  Ramberg-Osgood  parameters:  E = 10  psi, 

3 3 

Oq  7=24x10  psi,  n=6.66,  a^=16xl0  psi  and  v=0.33.  Present 
solutions  and  those  of  Refs.  [19,22]  are  shown  in  Figs.  8a 
and  8b.  The  finite  element  method  and  kinematic  hardening 
rule  are  used  in  Ref.  [19].  On  the  other  hand,  a finite  ele- 
ment formulation  together  with  an  isotropic  hardening  rule  are 
employed  in  Ref.  [22].  Good  agreement  is  observed  between 
present  solutions  and  those  of  [22].  At  higher  load  levels, 
results  of  Ref.  [19]  drift  away  from  present  solutions  and 
those  of  Ref.  [22] . 

Next  we  consider  a simply-supported  isotropic  spherical 
cap  under  a uniform  external  pressure.  Svalbonas  and  Levine 
predicted  an  elastic-plastic  buckling  load  for  such  a cap  in 
Ref.  [20].  The  material  under  consideration  obeys  a linear 
hardening  law  with  a ratio  of  tangent  modulus  to  Young's  mod- 
ulus of  0.1  and  yield  stress  to  Young's  modulus  of  0.002.  A 
shell  geometry  parameter  \=4  is  chosen.  Other  geometric  quan- 
tities are  given  in  Fig.  9.  Also  shown  in  this  figure  are 
the  present  load  vs  central  deflection  curve  and  that  of  [20], 
Comparison  between  these  two  sets  of  solutions  is  quite  good; 
the  present  analysis  predicts  a buckling  load  of  2050  psi 
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which  is  about  5%  below  the  value  of  2150  psi  estimated  in 

[20]  . 

The  final  example  problem  considered  is  a uniformly  load- 
ed clamped  spherical  cap  of  linear  hardening  material  (Fig.  10). 
Geometric  dimensions  and  material  properties  are  the  same  as  those  of  Fig.  9 
except  the  thickness.  Two  geometric  parameters  of  A=4  and  5.5  obtained  by 
varying  the  shell  thickness  are  analyzed.  Load  vs  central  deflection 
curves  are  plotted  in  Fig.  10.  Also  shown  in  this  figure  are 
the  results  obtained  in  Ref.  [6],  Good  agreement  between 
these  two  solutions  are  noted  for  the  spherical  caps  of  both 
A=4  and  5.5.  The  present  analysis  predicts  buckling  occurs 
at  load  of  2060  psi  for  A=4  and  1060  psi  for  A=5.5. 


> 
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ELASTIC-PLASTIC  BUCKLING  OF  SPHERICAL  CAPS  WITH  INITIAL 
IMPERFECTIONS 

Having  established  the  validity  of  plasticity  theory 
adopted  in  this  work,  we  turn  our  attention  to  the  main  ob- 
jective of  this  paper:  the  large  deformation  elastic-plastic 
buckling  analysis  of  spherical  caps  with  initial  imperfec- 
tions . 

Two  types  of  spherical  caps  (called  shells  A and  B in 
Fig.  11)  are  chosen  with  considerable  difference  in  geometrical 
dimensions.  Shell  A is  the  same  as  that  shown  in  Fig.  10; 
shell  B was  used  in  Ref.  [23]  for  the  dynamic  buckling  analy- 
sis. Three  geometric  parameters,  A=5,7.5  and  10,  for  each 
type  of  the  spherical  caps  are  selected  for  this  study,  the 
value  of  A is  obtained  by  varying  the  shell  thickness  h(A= 

2 [3  (1-v^)  ] (H/h)  . Material  properties  of  these  two 

spherical  caps  are  given  in  Fig.  11;  in  both  cases,  only 
linear  hardening  is  considered. 

The  axisymmetric  initial  imperfection  adopted  in  this 
study  is  of  the  dimple  type  which  was  also  used  in  Ref.  [8] . 
The  type  of  imperfection  is  expressed  mathematically  as 

2 ^ 

Wi  = (WiQ/b) (1-x  ) 

where  W^g  is  the  maximum  imperfection  which  occurs  at  the 
shell  apex.  Selection  of  this  expression  is,  in  fact,  quite 
arbitrary.  However,  it  does  provide  an  adequate  description 
for  actual  shells  since  the  important  parameter  is  the  maxi- 
mum eccentricity  and  not  the  imperfection  shape  function. 

In  the  buckling  analysis  of  spherical  caps  with  the  non- 
linearity in  both  geometry  and  material,  time  consuming 
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iteration  schemes  are  usually  involved,  and  hence  an  efficient 
numerical  procedure  is  highly  desired.  The  computer  program 
developed  in  this  study  provides  a restarting  capability  at 
any  load  level.  Thus,  in  the  buckling  load  calculation,  it 
is  able  to  halt  the  execution  at  any  load  level  and  restart 
from  that  point  by  a proper  choice  of  load  increments.  Also, 
it  is  able  to  restart  from  the  last  converged  load  level  by 
selecting  a smaller  load  step  after  the  iteration  fails  be- 
cause of  load  level  already  advancing  into  the  buckling  zone. 
With  this  facility,  more  accurate  buckling  loads  may  be  obtain- 
ed, and  considerable  computing  time  may  also  be  saved. 

A geometric  parameter  of  A=5  is  first  studied.  This  cor- 
responds to  h=0.26  in.  and  0.0227  in.  for  shells  A and  B, 
respectively.  The  magnitude  of  imperfections  considered  are 
W^g/h=0,  0.1,  0.5  and  1.  Results  associated  with  this  geo- 
metric parameter  are  depicted  in  Fig.  11  for  buckling  load  vs 
imperfection.  They  include  the  elastic  solution  and  those 
from  elastic-plastic  (linear  hardening)  analyses.  It  is  im- 
portant to  point  out  that  both  shells  A and  B yield  the  same 
elastic  solution.  However,  their  elastic-plastic  results  are 
different. 

For  the  purpose  of  references,  load  vs  central  deflection 
curves  are  plotted  in  Fig.  12  for  >=5  of  shell  A for  different 
imperfections.  Both  elastic,  elastic-plastic  solutions  are 
shown.  A redistribution  of  displacements  as  a result  of  the 
plastic  deformation  is  obvious  from  this  figure. 

Similar  to  what  is  shown  in  Fig.  11,  the  results  for  A= 

7.5  and  10  are  given  in  Figs.  13  and  14,  respectively.  Again, 
both  shells  A and  B provide  the  same  elastic  solution  for  each 
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of  A=7.5  and  10.  It  is  also  noted  that  both  solutions  from 
elastic  and  elastic-plastic  analyses  are  the  same  for  shell 
B in  Fig.  14.  This  means  that  the  shell  with  this  geometric 
parameter  has  the  thickness  so  thin  that  it  never  does  allow 
plastic  deformation  to  develop  before  buckling  occurs.  The 
effect  of  the  thickness  dimension  on  plastic  yielding  can  also 
be  seen  from  a comparison  of  results  among  Figs.  11,  13  and 
14:  the  larger  the  thickness,  the  more  the  influence  of  the 

plastic  deformation.  A summary  of  results  presented  in  Figs. 
11,  13  and  14  is  given  in  Table  1. 
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DISCUSSION  AND  CONCLUSIONS 

An  incremental  flow  theory  based  on  von  Mises  yield  sur- 
face and  in  conjunction  with  the  Prager-Ziegler  kinematic 
hardening  rule  is  utilized  for  the  calculation  of  plastic  de- 
formations. The  theory  is  then  incorporated  into  an  existing 
program  for  the  large  deformation  analysis  of  spherical  caps 
with  initial  imperfections.  The  problem  formulation  is  based 
on  equilibrium  equations  in  which  the  plastic  deformation  is 
treated  as  an  effective  plastic  load  and  lumped  together  with 
the  actual  applied  loads. 

In  every  solution  cycle,  the  material  property  and  effec- 
tive plastic  loads  are  held  constant,  and  the  problem  is 
reduced  to  an  elastic  one  whose  solutions  are  obtained  by  fi- 
nite differences  and  the  nonlinear  relaxation  technique  [17] . 
Effective  plastic  loads  are  then  modified  and  the  procedure 
is  repeated  until  they  correspond  to  the  computed  state  of 
stress  and  to  specified  stress-strain  relations  at  all  points 
over  the  shell  surface  and  through  the  thickness.  The  solu- 
tion procedure  mentioned  here  is  sketched  in  Fig.  5 for 
reference . 

It  is  noted  that  the  kinematic  hardening  rule  utilized 
in  this  work  accounts  for  the  Bauschinger  effect,  and  that 
perfectly  plastic,  linear  and  nonlinear  hardening  behavior 
are  all  included  in  the  program.  Nine  thickness  points  and 
a Simpson  rule  are  used  for  the  calculation  of  the  plastic 
deformation.  A dimple  type  of  imperfection  is  adopted  which 
provides  a quite  adequate  description  of  the  local  nature  of 
spherical  shells. 


In  general,  the  convergence  related  to  material  property 
iterations  is  quite  rapid.  For  example,  for  a load  level  not 
too  close  to  the  buckling  zone  and  with  a proper  load  incre- 
ment, only  about  2-5  iterations  are  usually  required  to 
satisfy  the  material  property. 

To  test  the  validity  of  the  plasticity  theory  implemented 
in  the  existing  program  [12] , five  elastic-plastic  plate  and 
spherical  cap  problems  were  solved.  Good  agreement  with  other 
solutions  was  obtained. 

Large  deformation,  elastic-plastic  buckling  loads  are 
calculated  for  two  types  of  spherical  caps  (Fig.  11) . For 
each  of  these  caps,  three  geometric  parameters  (A)  are  studied. 
Results  obtained  for  each  of  these  three  parameters  are  pre- 
sented in  Figs.  11,  13  and  14,  respectively.  From  these  re- 


sults, it  is  found  that  both  initial  imperfection  and  plastic 
deformation  have  the  same  effect  of  reducing  buckling  capa- 
city, and  that  as  the  magnitude  of  imperfection  increases,  the 
influence  of  the  plastic  deformation  becomes  less  significant. 
For  example,  at  W^g/h=l,  elastic  and  elastic-plastic  solu- 
tions are  almost  the  same  in  Figs.  11,  13  and  14.  It  is  also 
found  that  the  geometric  parameter  A,  which  is  used  as  an  im- 
portant factor  in  elastic  buckling  load  calculations,  becomes 
meaningless  for  the  elastic-plastic  buckling  analysis  of  these 
caps.  Elastic-plastic  buckling  loads,  in  fact,  depend  on  ac- 
tural  cap  geometric  dimensions  and  material  properties. 

As  a final  remark,  the  plasticity  theory  adopted  in  this 
work,  due  to  its  simplicity  in  the  mathematical  representation, 
is  quite  easy  to  apply.  More  importantly,  because  of  its  ca- 
pacity of  predicting  an  ideal  Bauschinger  effect,  the  theory 
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also  provides  very  realistic  results.  The  later  feature  is 
essential  for  the  analysis  of  cases  involving  reversed  or 
cyclic  loading  situations. 
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Fig.  2-  Flow  rule  and  kinematic  hardening. 


Fig.  3 - Kinematic  hardening  in  plane  stress. 
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Fig.  4- Nonlinear  stress-strain  curves 
(three  types  of  hardening). 
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Fig.  5 - Iteration  procedure  for  large  deformation  elastic - 
plastic  problems. 


r/a 

Fig.  6a-Deflection  profile  for  simply  supported  circular  plate  (Elastic 
ideally  plastic). 


Fig.  7a- Deflection  profile  for  uniformly  loaded  simply  supported 
circular  plate  (Elastic-linear  strain  hardening). 
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Fig.  7b - Distribution  of  circumferential  moment  for  uniformly  loaded 
simply  supported  circular  plate  (Elastic- 1 inear  strain 
hardening) 
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Fig.  8a  - Deflection  profile  for  uniformly  loaded  clamped  circular  plate 
(Elastic- nonlinear  hardening). 
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Fig.  8b-  Distribution  of  circumferential  moment  for  uniformly  loaded 
clamped  circular  plate  (Elastic- nonlinear  hardening). 
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Initial  Imperfection,  W^h 
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